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SUMMARY 


An improved photoelastic method has been developed for determining 
stress concentrations which occur in nonisotropic materials under load. 
Drucker*s method has been extended by means of an additional oblique 
photograph, and simpler equations for relating fringe order to principal 
stresses were obtained. The present method is considered more accurate 
than Drucker*s and also provides a means of determining the directions 
of the principal stresses. The determination of these directions, how- 
ever, is not sufficiently accurate, and the use of a stress coat was 
found desirable for greater accuracy. 


INTRODUCTION 


One of the pertinent problems involved in the use of nonmetallic 
materials for aircraft construction is the determination of the stress 
concentrations which occur in such nonisotropic materials under load. 

In an investigation of this problem the first type of stress distribu- 
tion selected for study was that which occurs around a loaded bolt 
passing through nonisotropic materials. In order to attack this problem, 
a photoelastic method has been developed. The present report describes 
this method and also its use in determining stress distributions in a 
loaded bolt. A review of the literature at the outset of the investi- 
gation indicated that the existing methods of photoelastic stress analysis 
were not entirely suitable for use in the problem under study. The usual 
and more simple photoelastic methods require the model to have a free 
edge, a requirement difficult to fulfill in this case. 

In order to overcome this limitation, several methods for analysis 
at single points have been developed (references 1 and 2), but these 
become tedious when analysis over an area is attempted. Others 
(references 3 and k) that are more suitable for analysis over an area 
require special apparatus as well as extreme care in the preparation 
of models. 

Drunker* s oblique-incidence method (reference 5) is essentially a 
point-by-point procedure unless his more complicated general equations 
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can be made applicable. Frocht states their use is practicable (refer- 
ence 6), but according to Drucker they cannot be applied accurately. 

It was found in the present study that Drucker* s method can be 
extended end made more useful by means of an additional oblique photo- 
graph and that simpler equations for relating fringe order to principal 
stresses can be obtained. The resulting procedure developed in the 
present investigation is considered more accurate than Drucker* s original 
one, and, in addition, makes available a means of determining the direc- 
tions of the principal stresses. The determination of these directions, 
however, is not sufficiently accurate, and the use of a stress coat 
(reference 7) was found desirable for greater accuracy. 

This work was conducted at the Forest Products Laboratory under the 
sponsorship and with the financial assistance of the National Advisory 
Committee for Aeronautics. 


APPAFATUS 


An apparatus was designed and built to utilize the double— oblique 
method of photoelastic analysis. The light, lenses, and camera for the 
polariscope were mounted on a track supported by a frame and bars (fig. 1). 
Each of the two side members of the frame was mounted on crossed pins 
that permitted the light to come from any direction within 0° to h5° of 
the normal position. 

The polariscope and supporting rig in position in a testing machine 
are shown in figure 1. The essential parts of the polariscope from left 
to right are transformer; mercury vapor lamp, filter, and 10— inch lens, 
the last three items covered by the shield; inch— diameter polarizing 

plate; quarter-wave plate; space for the specimen; quarter-wave plate; 
analyzing plate; 10-inch lens; and camera with 8- by 10— inch ground-glass 
viewing plate. These items are mounted on tracks that are pin— connected 
to the cross bars. The cross bar’s are free to pivot on the bolt hangers. 

A graduated • quadrant mounted, on the left end of the second cross bar is 
read against a pointer attached to the bolt hanger to determine the angle 
between the line through the horizontal pivots and the optic axis. The 
bolt hangers are suspended from channel iron members. Each channel is 
supported by a pin parallel to the hangers, and these pins are supported 
by alined horizontal pins. The angle between the horizontal plane and 
the projection of the optic axis on a plane perpendicular to the hori- 
zontally alined pins is read by means of a graduated quadrant mounted on 
the base and a pointer mounted on the horizontal pin. The supporting 
structure is made of pipe, channels, and bars welded together to leave 
access to the base of the testing machine and to the space for the 
specimen. Counterweights are placed where necessary to balance the 
assembly about the crossed pivots. 
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THEORY OF IMPROVED PHOTOEIASTIC METHOD 


Hp 




Drunker* o general equations are: 

-»n (1 + co3 2 0) cos 29 ± 2 cog 0 f n Q 2 - n n 2 sin 2 26j 


% 


( 1 ) 


2 sin 2 0 2 

-r^ (1 + cos 2 0) cos 26 t 2 cos 0 (n 2 - n 2 sin 2 n, 

= i L £ (2) 

2 sin 2 0 2 


where the symbols are defined as follows: 

np fringe order for the algebraically larger of the principal stresses 
fringe order for the other principal stress q 

fringe order from photograph with light normal to surface of model 
fringe order from photograph with light oblique to surface of model 


n q 

% 

n„ 


0 


angle between direction of light inside model and a normal to its 
surfaces 

orientation of p measured counterclockwise from x— axis 
Let % = np + n^; then, adding equations (l) and (2) gives 


n x “ 


(l + cos 2 0) cos 28 ±, 2 cos 0 ^n Q 2 ~ r^ 2 sin 2 26^^ 2 


sin 2 0 


(3) 


This equation may be put in the form: 


n^. 2 sin^ 0 + 2 x^ 1 ^ (l + cos 2 0) sin 2 0 cos 20 + r^ 2 sin* 0 cos 2 20 

+ ^(nn 2 - n-L 2 ) cos 2 0 = 0 (4) 


where n^ is the fringe number in the oblique photograph. 

An additional oblique photograph is taken with the light beam rotated 
so that the plane defined by its axis and the normal is at angle of 90° 
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to the similar plane associated with the first oblique photograph. Equa- 
tion (4) applies also to this photograph if 0 is replaced hy 0 + 

that is, the value of cos 20 is reversed in sign. On inserting n 2 , the 
fringe number from the second oblique photograph, in place of n^ and 
changing the sign of cos 20 , equation ( 4 ) becomes 


it* 2 s±n h 0 - 2n x n n (l + cos 2 0) sin 2 0 cos 20 + n, n 2 sir i k 0 cos 2 20 



( 5 ) 


Subtracting equation ( 5 ) from equation ( 4 ) leaves 

bn^ (1 + cos 2 0) sin 2 0 cos 20 + 4^n 2 2 - r^ 2 ) cos 2 0 = 0 

which can be put in the form: 

n x n^ 2 — cos 2 0 

Hn r^ 2 (l + cos 2 0) sin 2 0 cos 20 


( 6 ) 


( 7 ) 


The values of n^ r ip, and % at any point in an area may be 

determined from one group of photographs. The angle 0 is the angle 
between the oblique ray of light inside the model and a normal to the 
surface 3 of the model. The angle 0 is measured by using the stress- 
coat technique. The value of is determined by using these values 

in equation (7). It may be noted that no serious confusion of sign 
results from the use of this equation. The principal stresses are then 
easily determined: 


** “ Hp + n q 
“n - Hp - n q 


( 8 ) 


Addition gives twice the fringe number for one principal stress, and 
subtraction gives twice the fringe number for the other. Equation (7) 

becomes indeterminate when 0 = ^ and n-^ = rip, out this difficulty 
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may "be avoided "by talcing an additional pair of oblique photographs at an 
angle of 45° to the first pair. 

Equations (7) and (8) have the following advantages over equations (3) 
to (6) of Drucker*s paper (reference 5): 

1. Equation (7) is general; it applies for any position of the 
principal axes of stress. 

2. Drunker defines rip as algebraically larger than 
leaves some doubt regarding the identification of np and 
tions (7) and (8) make this identification clear. 

3. The right-hand member of equation (7) is independent of the sign 
of iijj, n-p and n^; thus the ratio n^/r^ has a unique value. This 

eliminates the confusion in the use of Drue ken* s equations because the 
sign of xIq is not known. 

4. Both methods are concerned with differences of fringe values 
measured from two photographs. The differences encountered by Drucker 
are iIq — 1^ cos 9 and n Q cos 9 — (in DruckBr*s notation). The 

differences encountered in the method described herein are n ^2 _ n^2 or 

( n l + n 2 ) (^ n l — Usually the value of . nq — ng is about double that 

of Hq — rijj f and the accuracy of the present method is therefore about 
twice that of .Drucker* s. 


n^. This 
n^ . Equa— 


COMPARISON OF RESULTS FROM PHOTQEIASTIC METHOD AND ELASTIC THEORY 


A 2-inch-diameter circular disk 0.252 inch thick was loaded in com- 
pression at opposite ends of a diameter. This model was placed so that 
its central plane was vertical and contained the line between the hori- 
zontal pivots of the polariscope Jig. Before load was applied, two 
photographs were taken with the light normal; one photograph is shown 
in figure 2 with white background and one in figure 3 with dark back- 
ground to show the initial fringe pattern. Load was then applied through 
round steel rods 0.125 inch in diameter. Photographs were taken with 
both white and dark backgrounds and with light direction (l) normal, 

(2) rotated vertically 45°, (3) rotated horizontally 45°, and (4) 
rotated 48° 19* from the normal in a 45° plane. These are shown re- 
spectively in figures 4 to 11. The photographs shown in figures 10 
and 11 were taken to provide data near the ends of the diagonal diamater 
of the model because of the possibility that data from the other photo- 
graphs might lead to indeterminate values of equation (7). This difficulty 
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vas not encountered in the analysis, and therefore the data from these 
photographs were not employed. After the fringe patterns under load 
had teen photographed, the load was removed and the photograph in 
figure 12 was taken to show the residual fringe pattern in the vicinity 
of the concentrated loads. 

The model was later sprayed with a stress— indicating lacquer. After 
the lacquer dried, load was again applied. Cracks perpendicular to the 
principal tensile stress formed in the lacquer. These were dye-etched 
to preserve them. After the coated specimen had remained under load 
about 1 hour, the lacquer and the model were chilled with air passed over 
dry ice. The load was immediately released; the result was a formation 
of cracks perpendicular to the principal compressive stress. These also 
were dye-etched. A photograph of the etched cracks is shown in figure 13. 

A circle 1.20 inches in diameter had been drawn on each surface of 
the disk so that correct measurements could be made on the prints in 
units of the radius of the scribed circle. A transparent interpolator 
was used to read the location of fringes to a radius of 0.01. Curves 
with fringe number as ordinate and distance from the center of the model 
as abscissa were prepared for the vertical diameter of the disk (fig. 1^) 
the horizontal diameter (fig. 15), and a diameter at h5° to the vertical 
diameter (fig. 16 ) . The values of the angle 9 between direction of 
principal stress at various points along a h5° diameter and the x— axis 
were determined from the photograph in figure 13 and are plotted in 
figure IT. The index of refraction was obtained from the rotated views 
showing the two scribed circles displaced in the horizontal and in the 
vertical directions. 

From the values of n^, n n , 6, and 0 so obtained, the 

values of n^ were computed (by equation (7)) for several points on 
the vertical, horizontal, and 1+5° diameters. From n^ and n^, the 
values of np and n^ were obtained and are plotted in figures 18 and 19 

for the horizontal and vertical diameters, respectively. A smooth curve 
was drawn through these values in each figure. 

The mathematical determination of the principal stresses is based on 
the theory given by Timoshenko (reference 8) from which the following 
general equations for the principal stresses at any point in a disk under 
edgewise compression with point load were derived: 


°u = a - 



■) 


( 9 ) 


cr T = a 


( 10 ) 
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Factors in equations ( 9 ) and (10) are: 

Yg “ |( x / d ) 2 ( 2 - cos 2A) + (x/d)^ 


cos 2 B 


^ — ^(x/d ) 2 cos 2A + (x/d)^ 


(ID 


cTt = a 


| - I coo A 


(x/d ) 2 — ^ cos A + 1/4 
d 


( 12 ) 


cr 2 = a 


lx . 
■k + — cos A 

v Cl 


(x/d ) 2 + i cos A + 1/4 
& 


(13) 


where 


V a v principal stresses at a point 
x radial distance from center of disk to point 

d diameter of disk 

A angle between loaded diameter and line from center of disk to 

point 



P load 


From equations ( 9 ) and (10) the following formulas were obtained for 
the principal stresses on the horizontal and vertical diameters: On the 

horizontal diameter with 0 = £ and cos 0=0, 


u 


= a 



(d 2 + 4x 2 ) 2 


= a 


1 


d 2 x 2 




(14) 


(15) 
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On the vertical diameter with 0=0 and cos 0=1, 



The numerical value of a required to express the theoretical stresses 
in units of fringe order was determined by the method of least squares by 
using the crosses on the curves in figures 18 and 19 for the horizontal 
and vertical diameters, respectively. This value of a was used to 
determine the principal stresses in units of fringe order on the diameter 
at U-5° to the horizontal, on the horizontal diameter, and on the vertical 
diameter. Points were not taken at the center of the disk because the 
wide fringes were difficult to locate accurately; points were not taken 
near the edges to avoid the initial fringes (figs. 2 and 3 ) and. the 
residual fringes (fig. 12). 

In figure 20 the computed principal stresses in fringe-order units 
are shown by the curves, snd the experimentally determined principal 
stresses in the same units are shown by the plotted points. 


DISCUSSION 


A comparison of the stress distribution obtained by the photoelastic 
method with that obtained by the mathematical theory of elasticity is 
of interest. Such a comparison is given in figure 20, in which curves 
obtained by the theory of elasticity and points obtained by the photo- 
elastic method are compared. The fringe value, which relates the stress 
to the fringe order, was obtained by fitting the mathematical curves to 
the experimental data of figure 20(a) to 20(d) by the method of least 
squares as previously described. The fringe value so determined was used 
in drawing the curves in figure 20(e). 

The data agree with the mathematical curves reasonably well in most 
instances. The deviations can readily be explained by an analysis of 
the experimental technique. Of course, it was impossible to apply true— 
point loads to the disk. Also the stresses near the points of applica- 
tion wei-e of necessity extremely high, and the material of the model 
became plastic at these points. The result was the development of local 
areas of approximate hydrostatic stress which increased the horizontal 
tension in adjacent parts of the specimen. The pr®ence of such tensile 
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stresses is indicated in figure 20(a). These plastic areas also tend to 
decrease the vertical compressive stress at the center of the specimen 
and to increase it near the ends of the horizontal diameter. (See 
fig. 20(d).) Of course, the photoelastic method is not expected to 
yield accurate stresses at locations in the speciiren at which plastic 
flow is present; thus the measured stresses at or near the residual 
fringes shown in figure 12 are in error. This error is shown in 
figure 20(c) and 20(e). 

The data taken on the 45° diagonal of the disk serve "best as a 
verification in the improved photoelastic method, since it was necessary 
to use the isostatic lines obtained by the stress-coat method (fig. 13 ) 
to determine the experimental values. Although the principal stresses 
on the vertical and horizontal diameters could have been obtained by 
Drucker*s method, the improved method was used, with 0=0 or 

because of the simplicity of its application. On the 45° diameter, the 
ratio of the photoelastically observed principal tensile stress to the 
computed value is 1.00 at the center of the disk and 1.10 at the position 0.4 
of the radius from the center. This ratio for the compressive principal 
stress is 0.9 6 at the center of the disk and 0.94 at 0.4 of the radius 
from the center. Stresses closer to the edge of the disk are considered to 
be inaccurate because of the errors previously discussed. 


CONCLUDING REMARK 


It can be concluded that the improved oblique method of photoelastic 
analysis yields reasonably accurate results over regions of the photo- 
elastic model which are not subjected to excessive stress concentrations. 


Forest Products Laboratory, Forest Service 
U. S . Department of Agriculture 
Madison, Wis., July 15, 19^6 
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Figure 2.- Disk before load was applied. Light normal, white background. 
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Figure 4.- Disk under load. Light normal, white background. 
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Figure 5.- Disk under load. 


Light normal, dark background. 
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Figure 6.- Disk under load. Light deflected 45°, white background. 
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Figure 7.- Disk under load. Light deflected 45°, dark background. 
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Figure 8.- Disk under load. Light rotated 45°, white background. 
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Figure 9.- Disk under load. Light rotated 45°, dark background. 
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Figure 10.- Disk under load. Light oblique 48° 19', white background. 
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Figure 11.- Disk under load. 


Light oblique 48° 19', dark background. 
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Figure 12.- Disk after load was removed. 


Light normal, white background. 
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Figure 13.- Etched cracks showing directions of stress in lacquer applied 

to face of disk. 
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AVERAGE DISTANCE FROM CENTER , RADII OF SCRIBED CIRCLE 

Figure 14.- Fringe -order curve for vertical diameter of disk. 
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Figure 15.- Fringe-order curve for horizontal diameter of disk. 
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AVERAGE DISTANCE FROM CENTER , RADII OF SCRIBED CIRCLE 

Figure 16.- Fringe-order curve for diameter at 45° to vertical diameter of disk. 
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PRINCIPAL STRESSES , FRINGE- ORDER UNITS 



Figure 18.- Principal stresses on horizontal diameter. 
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Figure 19.- Principal stresses on vertical diameter. 
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Figure 20.- Comparison between computed and experimentally 
determined principal stresses. 


